cpi uiaii 



Euclidean matrix theory of random lasing in a cloud of cold atoms 



A. Goetschy 1 and S.E. Skipetrov 1 



O 

> 
O 

(N 
C 

• i-H 

o3 



i 

C 

o 
o 



> 



(N 
O 



X 



Universite Grenoble 1/CNRS, LPMMC UMR 5493, B.P. 166, 38042 Grenoble, France 
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Abstract. - We develop an ab initio analytic theory of random lasing in an ensemble of atoms that 
both scatter and amplify light. The theory applies all the way from low to high density of atoms. 
The properties of the random laser are controlled by an Euclidean matrix with elements equal 
to the Green's function of the Helmholtz equation between pairs of atoms in the system. Lasing 
threshold and the intensity of laser emission are calculated in the semiclassical approximation. 
The results are compared to the outcome of the diffusion theory of random lasing. 



Introduction. — We call 'random' a laser in which 
the feedback is provided by the multiple scattering of light 
on the random heterogeneities of the active medium and 
not by a well-defined external cavity pQ. Recent theo- 
retical models of random lasers rely on expansions of the 
laser field in terms of overlapping modes of 'random cav- 
ities' formed by these heterogeneities [2J[3]. Alternative 
approaches consist in solving Maxwell-Bloch equations nu- 
merically [IHS] or within the diffusion approximation [THS] ■ 
The latter has the advantage of yielding a simple criterion 
for the lasing threshold but it lacks rigorous justification 
and breaks down in the strong scattering regime pQ. 

In the present Letter we develop a new approach to the 
problem of random lasing that does rely neither on the ex- 
pansion of the laser field in terms of cavity modes, nor on 
the diffusion approximation. It is based on our recent ana- 
lytic results for the 'random Green's matrix' that describes 
propagation of light between scattering centers in a ran- 
dom medium [10] . To demonstrate the power of this new 
approach, we consider random lasing in an ensemble of a 
large number N of identical atoms in free space, a prob- 
lem of recent interest [9]|TT]|T2] . Dynamic equations for N 
atoms that both scatter and amplify light are presented 
and analyzed in the semiclassical limit. We obtain analytic 
results for the lasing threshold and the average emitted in- 
tensity, thus achieving an important progress with respect 
to previous works on similar systems by Savels et al. (who 
treated lasing in ensembles of N < 5 three-level atoms) 
[12] and Froufe-Perez et al. (who dealt with N 3> 1 two- 
level atoms but in the diffusion approximation) [SJ. Our 
approach can be extended to deal with more 'standard' 
random lasers in which scattering centers ('particles') are 



embedded in an amplifying homogeneous matrix [14] . 

The model. — Consider a gas of N three-level atoms 
at random positions r, (i = 1, ...,N) in free three- 
dimensional space. The atoms are subject to a strong 
external pump field of amplitude E p (r), resonant with 
the transition from the ground state \ai) to the upper 
auxiliary level |cj) of each atom. They then rapidly de- 
cay to the upper level |6j) of the laser transition at a 
rate r c f, ^> Fb a = F ^ F ca . Interaction of atoms with 
the electromagnetic field which is near-resonant with the 
transition from \bi) to |o») (energy difference fiwo) can be 
described using standard approaches [15], although cer- 
tain subtleties should be treated with care. Relegating 
the derivation of dynamic equations for population imbal- 
ances Hi — \bi}(bi\ — \ai}(ai\ and atomic lowering operators 
Sj~ = |di)(6i| to a future publication [14] , we present here 
their semiclassical approximation: 

^ = _(i + w i )n i + wi-i + i[s?n i -c.c.], (ia) 



dSj 
dt 



Si - -HA, (lb) 



where time t is in units of T -1 . IF = (0|II,*|0) and Si = 
(OlS^lO), with |0) the vacuum field state [13], measure the 
difference between populations of levels \bi) and \ai) and 
the (induced) dipole moment of the atom i, respectively. 
The pumping rate is given by Wi = [dj • E p (ri)] 2 /H 2 T c bT 
with di the dipole moment for the \ci) — > \ai) transition, 
fii is the i-th element of the vector Q = G(wo)S, with 
S = (Si, . . . , Sn)- It is proportional to the amplitude of 
the electric field at atom i and is due to radiation of all 
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other atoms. The N x N non-Hermitian Green's matrix 
G with elements 

G.M = (l-%)g^ (2) 

couples different atoms (fco = wo/ c = 27r/Ao). This matrix 
describes propagation of light between atoms and belongs 
to the family of Euclidean random matrices defined by 
positions of N points in the Euclidean space jTB] . 

Equations (JTJ) are derived in the scalar approximation 
for the electromagnetic field and assuming V <C wq, c/R 
(with R the size of the atomic cloud). They can be 
regarded as a generalization of the optical Bloch equa- 
tion |15| to an ensemble of identical, incoherently pumped 
atoms. In the absence of coupling between atoms, they de- 
scribe an isolated atom and have the stationary solution 
11, = (Wj — l)/(Wi+l), Si = 0, that shows that population 
inversion IL, > can be achieved for Wi > 1. This thresh- 
old for achieving population inversion, as well as the power 
broadening of the transition (the natural linewidth T is in- 
creased by a factor 1 + Wi ) , are due to sharing of the same 
ground state by the pump and the lasing transition. The 
dimensionless polarizability a = (k^/AiT)a = Si(ui) /fij(u;) 
of an atom pumped at a rate W follows from eqs. ([T]) (see 
also ref. [12]): 

W-l 172 
& ^ W) = WTT (u-. Q ) + l (l + W)T,2 - (3) 

Lasing threshold. — When coupling between differ- 
ent atoms is at work, the stationary solution Si = of eqs. 
|T]) may lose its stability for a sufficiently strong pump. 
Following standard semiclassical laser theories [17] , we will 
associate this instability with reaching the lasing thresh- 
old. The stability analysis of nonlinear eqs. ([1]) shows that 
for uniform pump Wi = W lasing starts when at least one 
of the eigenvalues of G satisfies 

2W 

— ^ImAfe > (l + W0+ImA fc . (4) 

The left-hand side of this condition can be regarded as gain 
that depends on both the pumping rate W and scattering 
(through Afe), whereas the right-hand side contains pump- 
dependent losses due to spontaneous emission (1 + W) and 
leakage out of the system (ImAfe). As counterintuitive as 
it may seem, it follows from eq. ((4]) that random lasing 
takes place when ImA^ (that quantifies losses due to open 
boundaries in the absence of pump) exceeds (1+W) 2 / (W— 
1) and W >1. 

We now propose a less technical and physically transpar- 
ent derivation of the threshold condition ((4} which, in ad- 
dition, can be generalized to an arbitrary pumping scheme. 
On the one hand, the electric field at atomic positions is 
Sl(to) = G(uj)S(uj), whereas on the other hand, the linear 
response of atoms to the field implies S(w) = A(uj)fl(co), 
where A is the diagonal matrix A — diag[<Si] and on is 



the dimensionless polarizability of atom i. It then fol- 
lows immediately that the linear description breaks down 
and lasing starts when one of the eigenvalues of the ma- 
trix product G(u>)A(uj) is equal to one. For atoms in free 
space, G(uj) can be safely replaced by G(ujq). Therefore, 
for identical atoms and uniform pumping (<S,; = a) the 
threshold condition reduces to 

A fe = 4. (5) 
a 

This equation illustrates that laser threshold results from 
an interplay of single-atom properties (described by the 
polarizability a) and geometry-dependent collective effects 
(quantified by the random eigenvalues Ak of the Green's 
matrix G). If we substitute the polarizability ^ into eq. 
([5]), we recover eq. (@|. But eq. ([3]) is more general than eq. 
(HJ) and is not restricted to lasing in a system of three-level 
atoms. We can also apply it, for example, to an ensemble 
of two- level atoms (resonant frequency ojq) in the field of 
a strong near-resonant coherent pump (frequency ujq + A, 
Rabi frequency Q). When illuminated by a weak probe 
light at a frequency cj + A + S, each atom behaves as 
if it had an effective polarizability a(6, A, fi) with a long 
but manageable explicit expression that can be found in 
Refs. 9,18.19 . Optical gain in such a system is sometimes 
referred to as 'Mollow gain' [TS] . 

An easy way to visualize the threshold condition ([5]) is 
to draw the two-dimensional domain T>\ occupied by the 
eigenvalues of G and the region T> a spanned by 1/a when 
its free parameters — oj and W in the case of eq. ((3]) — are 
varied, on the complex plane. Random lasing takes place 
when 2?a and T> a touch (threshold) or overlap. This is 
illustrated in fig. [1] for N 3> 1 atoms in a sphere of radius 
R 3> Ao- In this figure, we adjusted the parameters for 
the random laser to be slightly above threshold: T>\ and 
T> a barely overlap. Whereas V a is easy to determine when 
a is known as a function of its parameters, finding 2? a is 
much less trivial. Here we make use of our recent results 
for the eigenvalue distribution of the Green's matrix @ in 
the limit of large N [10] . The distribution and the bound- 
ary of its support X>a on the complex plane depend on 
two dimensionless parameters: the number of atoms per 
wavelength cubed pAp 1 and the on-resonance optical thick- 
ness bo = 2R/£, where p is the number density of atoms 
and i = Kq/ Ak p is the on-resonance scattering mean free 
path in the absence of the pump. The eigenvalue domain 
X>a consists of two parts: a (roughly circular) 'bulk' and a 
pair of spiral branches [TU1I20] . Depending on the partic- 
ular model of atomic polarizability a, either the bulk or 
the branches may touch V a , as we now discuss. 

We first focus on the lasing threshold due to the bulk 
of eigenvalues. Combining the analytic equation for the 
borderline of T>\ at low density pA§ < 10 10J and eq. 
([5]) results in a threshold condition that depends on the 
optical thickness &o but not on the density pAg: 

h \&\ 2 h (h Im&) = 1, (6) 
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Fig. 1: The domain D a (hatched) spanned by \/& and the domain T>k (blue area delimited by the solid line) occupied by the 
eigenvalues A of the random Green's matrix ((2}. (a) Incoherent gain &(uj, W), see eq. {SJ. (b) Coherent Mollow gain a(S, A, Q) 
with A = r [Dl l 18lfT5] . Lasing occurs when T> a and 7?a overlap: regions (1), (2), (3). The borderline of T>a is given by eq. (JU| 
with the optical thickness bo = 40 in (a) and bo = 140 in (b). The dashed lines show the borderline of T>a following from the 
diffusion approximation [eq. (0]. 



where h(x) = [3 - 6a; 2 + 8x 3 - 3(1 + 2x)e- 2x ]/Qx A . Note 
that for both gain mechanisms considered in this Letter, 
the threshold condition ^ involves the eigenvalue with 
the largest imaginary part, as can be seen from fig. [T] 
We calculated (max(ImA)) based on our non-Hermitian 
random matrix theory jTU] and found excellent agreement 
with numerical results, see fig.^a). It is quite remarkable 
that the agreement is present at all values of parameters, 
including high densities pXq 1 that were necessary to 
reach large optical thicknesses bo 3> 1 in numerical calcula- 
tions with moderate N < 10 . Because it is (max(ImA)) 
that controls the laser threshold, we conclude that our 
theory applies to random lasing all the way from weak 
(p\y -C 1) to strong (pAg 3> 1) scattering regime. 

It is interesting to compare the threshold condition © 
with the one obtained in the diffusion approximation. The 
latter amounts to solve the diffusion equation for the av- 
erage intensity of light with gain included as a negative 
absorption [7HS] ■ The threshold is reached when the solu- 
tion diverges. This yields the following threshold condition 



V3 
2tt 



Ima 



1 



b \a 



~.\2 



1. (7) 



This equation is similar to our result ([6|) at large opti- 
cal thickness b — bo\a\ 2 3> 1 but deviates significantly at 
b < 1, as can be seen from fig.[TJ Consequently, the predic- 
tions of eq. (J7J for the laser threshold (that is reached at 
b < 1, see fig. HJ turn out to be inaccurate. In particular, 
our eq. ([5]) predicts that the minimum on-resonance op- 
tical thicknesses required for random lasing are &ocr — 35 



1 We use the value of the extrapolation length zq = 2/3 instead 
of z = 0.71 in Ref. [5]. 



for the incoherent and 6g cr — HO for the coherent pump. 
This is significantly less than 50 and 200, respectively, fol- 
lowing from eq. ((7} . Analysis of the right eigenvectors R*. 
(modes) of the matrix G shows that at all densities pAg, 
the mode that reaches the threshold first is extended over 
the whole atomic cloud [see fig. [3[c)], even when the sys- 
tem may support localized modes as well [see modes (b) 
and (d) in fig. [3J. This is specific for the models considered 
here in which, in particular, scattering and gain are due to 
the same atoms, and in contrast with systems where gain 
and scattering are independent and (pre-)localized modes 
may be better candidates for lasing |21[l22] (see ref. [6] for 
a recent review). 

In the high-density limit p\g — > oo, the eigenvalues of 
G that have large imaginary parts collapse on a line de- 
scribed by a simple equation [10) which, combined with eq. 
([5]), yields the lasing threshold condition for a continuous 
medium with a refractive index n(a) = (1 + &pXq/27t 2 ) 1 ^ 2 : 



n(a) — 1 



n(a) + 1 



3 4in(ci)fco^ 



1. 



(8) 



In this limit the problem looses its statistical nature and 
the random laser turns into a 'standard' laser with the 
feedback due to (partial) reflections at the boundaries of 
a homogeneous amplifying medium. 

Let us now analyze the role of the spiral branches of 
T>a- As we illustrate in fig.[3][b), the eigenvalues belonging 
to these spirals correspond to eigenvectors (modes) local- 
ized on pairs of very close points |r-j — r^| -C A [TU1I2T)] . 
These are the super- and sub-radiant states of a pair of 
atoms. For the uniform incoherent gain ©, the branches 
do not overlap with T> a [fig.QJa)], whereas the lower, 'sub- 
radiant' branch overlaps with T> a for the coherent Mollow 
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Fig. 2: Mean maximum value of the imaginary part (a) and mean minimum value of the real part (b) of eigenvalues A of 
the N x N random Green's matrix ([2]). Our analytic results (solid lines) are compared with numerical simulations for three 
different matrix sizes iV (symbols). The horizontal dotted lines indicate the absolute lasing thresholds: (1) for the incoherent 
gain [eq. (|3}] and (2) and (3) for the coherent Mollow gain [91 11811X9] . In (a), analytic results depend on both bo and pA(j 
[10], except for pA(j < 10 when they reduce to eq. (|6]) with a replaced by 1/i (max(ImA)) (dashed line). The analytic result 
(min(ReA)) = -r(2/3)(ATpA 3 ,/127r 2 ) 1/3 [solid line in (b)] is valid for pA 3 , < 10. 




Fig. 3: (a) Eigenvalues A of a single random realization of 
the Green's matrix G (dots) for a cloud of optical thickness 
bo = 40, composed of N = 10 3 atoms, (b)-(d) Intensities \R\\ 2 
corresponding to the mode in the subradiant branch, localized 
on a pair of atoms (b), the mode with the largest ImA (c) 
and the mode corresponding to the smallest |A| (d). A mode 
Rfe = {Rl , i? 2 . , . . . , Rf* } is represented by spheres centered at 
positions of atoms rj and having radii equal to lx (b), 100 x 
(c), and 10 x |i?^| 2 (d). 



gain [fig. [IJb), region (3)]. Thus, in the latter case the 
solution Si — of eqs. (QJ may lose its stability due to 
the eigenvalue with the smallest real part belonging to 



this branch. We calculated (min(ReA)) from the proba- 
bility to find two points separated by a distance Ar > xR: 
p= [l-x 3 (l~9x/l6 + x 3 /32)] N( - N - 1 ^ 2 . The result scales 
with (A^pAg) 1 / 3 , in good agreement with numerical sim- 
ulations as long as pXf, < 10 [fig. |2{b)]. When (— ReA) 
exceeds a critical value, a pair of closely located atoms on 
which the eigenvector (mode) associated with the eigen- 
value A is localized, starts to emit coherent light. On 
average, the threshold for this effect is given by the con- 
dition 

where T(x) is the Gamma function. Formally, this emis- 
sion of light by a pair of pumped atoms may be called 
'laser', especially given the fact that one-atom cavity lasers 
[23] and few-atom random lasers [13] were already dis- 
cussed in the literature. It is very different from the col- 
lective laser mechanism leading to cqs. ([5]) and © and as- 
sociated with eigenvectors extended over the whole atomic 
system. Whereas eqs. ^ and ([5]) are good estimates of 
the threshold even for a single atomic configuration, the 
threshold for light emission by a pair of atoms is expected 
to fluctuate strongly around its typical value given by eq. 
(|9]) . In addition, this emission is strongly affected by quan- 
tum effects not included in our analysis. In fact, quantum 
fluctuations wash out the sharp threshold obtained in the 
semiclassical framework (14] . 

Laser intensity. — Let us now study the dynamics of 
laser emission slightly above threshold. In the vicinity of 
threshold, IL can be adiabatically eliminated from eqs. ([TJ 
[17] . Keeping only the lowest-order nonlinear terms, we 
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Pump W Pump W 



Fig. 4: (a) The minimal (critical) optical thickness 6ocr necessary for lasing, following from our Euclidean matrix theory ((6| 
(solid line) and from the diffusion approximation {7J (dashed line), (b) The average stationary intensity (7) at bo = 40 obtained 
from the numerical solution of eqs. Q} for N = 10 3 after averaging over 10 random configurations of atoms (symbols). For each 
configuration and each W, we averaged the numerical solution I(t) over t = (200-250)r~ 1 . The analytic solution (j 12|) is shown 
by the solid line. The dashed line corresponds to yo = 9.15 instead of 8.93 for the solid line. 



obtain an equation for the dimensionless field amplitude 
fi = GS: 



~dt 



>.\Y + ga-gc\vl\ 



GBG- 1 



n, (io) 



where we introduced N x N diagonal matrices A = 
diag[(Wi - l)/(Wi + l)]/2, B = diag[Wi + l]/2, C = 
diag[(Wi - l)/(Wi + l) 3 ], and \(l\ 2 = diag[|^| 2 ]. Note 
that the eigenvectors of the linear kernel of eq. (|10[) co- 
incide with the right eigenvectors Rfc of the matrix G, 
that play the role of eigenmodes of the 'cold cavity', only 
if the pump is uniform: Wi = W. Restricting further 
consideration to the latter case, we introduce an ansatz 
ft = ^2 k a k {t)e~ lUkt Tl k [T7] and obtain rate equations for 
mode intensities I k — |afc| 2 (Rfc|Rfc): 



dh 
dt 



= -2n k I k + ^ Wnklnh 



(11) 



where 2n k = T[(W + 1) - ImA k (W - 1)/(W + 1)], W nk = 
-4rim(A fe 7„ fe )(W / - l)/(W + l) 3 , and a free-running sit- 
uation (no phase locking) is assumed. The threshold for 
the mode k is given by the condition n k = and depends 
only on the eigenvalue , whereas the mode competition 
above the threshold involves the overlap of eigenvectors 
lnk = £ti(^D 2 KlVEtiKI 2 - K ^ worth noting 
that although rate equations similar to eq. (fTTj) appeared 
in previous works on random lasers [2], loss rates K k and 
nonlinear couplings W n k were most often assumed to fol- 
low from ad hoc random matrix models, except in one- 
dimensional systems where they could be calculated with 
a reasonable effort. We, in contrast, provide explicit gen- 
eral expressions for these quantities and show that they are 
determined by the eigenvalues and eigenvectors R^ of 
the random Green's matrix. The link between n k , Wnk 



and Afc, Rfc is independent of the geometry or dimension- 
ality of the problem. 

At low atomic density pAg < 10, R\ behave almost as in- 
dependent Gaussian random variables and thus N("f nk ) ~ 
1 + 25 nk . This allows us to calculate the average intensity 



in the stationary regime (I) = (J2iLi 1^ 



(E fc 4): 



(I)=N 



{l + Wf 



1 



(W + l) 2 1 
W-l w 



(12) 



where y is the solution of y 2 = (3 /8)b h(—b /2y ) ■ This 
result is in good agreement with the numerical solution 
of eqs. ([IJ, as we show in fig. @|b). The average number 
of lasing modes, i.e. of eigenvectors Rfc of G that have 
non- vanishing amplitudes a k (t) in the expansion of O, in 
the long-time limit, is (N L ) = (2AT) 1 / 2 [(t^ - 1)/(W + 
I) 2 — l/yo] 1/ ' 2 2/o/(2/o + !)• Note that in the absence of 
mode competition, (I) and (Nl) would scale as A^ 2 and 
N, respectively. 

An interesting feature of lasing in a cloud of cold atoms 
illustrated by fig. HJb) is the halt of lasing at too strong 
pumps. This can be easily understood by noting that ran- 
dom lasing requires both amplification and scattering to 
be sufficiently strong, and that both of these important in- 
gredients are provided by the same atoms. At low pump 
[W < 2 in fig. Hfb)], the scattering is strong, but the 
amplification is not enough to lase. In contrast, when the 
pump is strong (W > 5), the scattering strength decreases 
because the atomic transition starts to be saturated, and 
lasing stops. 

Possible extensions. — The results presented above 
concern the simplest case of random lasing in an ensem- 
ble of atoms in free space, in the scalar approximation 
and under a spatially uniform pump. This allows us to 
demonstrate the power of the Euclidean matrix approach 
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which, in this case, yields explicit analytical results for 
the lasing threshold [eqs. ©, flSJ) and ©] and the laser 
intensity [eq. (|12l) ]. However, the approach is not lim- 
ited to this simplest case and can be extended to take 
into account additional complications that may arise in 
an experiment. If, for example, the pump is nonuni- 
form, finding the lasing threshold will require diagonal- 
ization of the matrix GA — iGBG^ 1 instead of G, as fol- 
lows from eq. (|10[) . For a complex distribution of pump, 
this can be realized only numerically. Another impor- 
tant situation is lasing in an ensemble of passive scat- 
terers embedded in an amplifying matrix with an am- 
plification coefficient /i. In this case, one has to put 
W = and study the eigenvalues of the Green's ma- 
trix G,j(w ) = (1 - <%) exp[(ifc + //)|r f - |]/fc |r-i - rj\, 
which can be done analytically [Tl]. Finally, taking into 
account the vector nature of light requires using the dyadic 
Green's function instead of the scalar one and working 
with 3N x 3A^ matrices. The threshold condition is still 
given by eq. ((5|) and at low atomic densities considered in 
this work, the analysis of the threshold due to the bulk of 
eigenvalues can be reduced to the scalar approximation. 
However, the spiral branches of the distribution (see fig. 
Q} are sensitive to the vector nature of light and thus the 
subradiant laser threshold (JH) will be modified. 

Conclusion. — Quite remarkably, the general frame- 
work developed in this paper and, in particular, the 
threshold condition ([5]), apply for any dimensionality of 
space, any atomic polarizability, any number and concen- 
tration of atoms, and any form of the Green's matrix G(ui) 
that, in particular, can account for an external cavity and 
amplification or absorption of light in the space between 
the atoms. This opens a way to treat various random 
laser systems analytically without strong approximations 
or heavy numerics. 

Application of the Euclidean random matrix theory to 
the problem of random lasing in a cloud of cold atoms 
allowed us to find the lasing threshold without relying on 
the diffusion approximation or transport theory, as well as 
to study the intensity of laser emission beyond threshold. 
We predict the possibility of random lasing in a cloud of 
cold atoms for on-resonance optical thickness exceeding 
35 for three- and 110 for two- level atoms ('Mollow laser'). 
We show that mode competition plays an important role 
in the random laser and leads to the scaling of the number 
of modes participating in lasing with \J~N (instead of N in 
the absence of mode competition) , where N is the number 
of atoms. At the same time, the laser intensity scales with 
N (instead of iV 2 ). 

Finally, our approach based on Euclidean random ma- 
trix theory may be also applied to understand the role of 
disorder in other collective effects that take place in en- 
sembles of cold atoms, such as the collective spontaneous 
emission |24) or the radiative force experienced by an op- 
tically thick atomic sample |25j . 
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